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We present the new paradigm of critical current by design. Analogous to materials by design, it 
aims at predicting the optimal defect landscape in a superconductor for targeted applications by 
elucidating the vortex dynamics responsible for the bulk critical current. To highlight this approach, 
we demonstrate the synergistic combination of critical current measurements on commercial high- 
temperature superconductors containing self-assembled and irradiation tailored correlated defects 
by using large-scale time-dependent Ginzburg-Landau simulations for vortex dynamics. 
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1. INTRODUCTION 

The interaction of vortex matter with defects in ap¬ 
plied superconductors directly determines their current 
carrying capacity. Defects range from chemically grown 
nanostructures and crystalline imperfections to the lay¬ 
ered structure of the material itself. The vortex-defect 
interactions are non-additive in general, leading to com¬ 
plex dynamic behavior that has proven difficult to cap¬ 
ture in analytical models. With recent progress in ad¬ 
vanced computing, a new paradigm has emerged that 
aims at simulation-assisted design of defect structures 
with predictable powers. This concept of critical current 
by design aims at predicting the optimal defect land¬ 
scape for targeted applications by elucidating the vor¬ 
tex dynamics responsible for the bulk critical current. 
To highlight this approach, we demonstrate the syner¬ 
gistic combination of critical current measurements on 
commercial high-temperature superconductors contain¬ 
ing self-assembled and irradiation tailored correlated de¬ 
fects with large-scale time-dependent Ginzburg-Landau 
numerical simulations of vortex dynamics. The qualita¬ 
tive agreement between experiments and simulations al¬ 
lows making numerical predictions for a wider parameter 
range and hence opens the route to the critical-current- 
by-design paradigm. 

Superconducting vortices appear inside superconduc¬ 
tors in the presence of sufficiently large magnetic fields 
and are responsible for the entire electromagnetic behav¬ 
ior of applied superconducting materials. In particular, 
the onset of vortex motion limits the highest dissipation¬ 


less electrical current — the critical current — the su¬ 
perconductor can carry. At the microscopic level, vor¬ 
tices can be viewed as non-superconducting cores that 
are surrounded by circulating persistent superconducting 
currents. The radius of the core is approximately given 
by the superconducting coherence length, which for 
high-temperature superconductors is of the order of 1.5- 
3 nm, and the circulating supercurrents extend out to the 
magnetic penetration depth. A, which is of the order of 
100-200 nm. When an external electric current is passed 
through the superconductor, the vortices experience the 
Lorentz force, which sets them in motion, resulting in a 
voltage drop across the superconductor and hence dis¬ 
sipation: the superconductor looses its superconducting 
properties.^ A variety of defects in the superconducting 
material (precipitates, point defects, grain boundaries, 
dislocations, stacking faults, strain fields, etc.) can in¬ 
duce spatial variations of the superconducting state that 
serve to pin the normal vortex cores and prevent their 
motion. Thus, an important goal of materials science 
of applied superconductors is the design and synthesis 
of the most effective vortex pinning microstructures that 
results in the largest critical current. Although a tremen¬ 
dous amount of knowledge, systematics and practical ex¬ 
perience has been gained in this pursuit, see for instance 
reviews,the fundamental solution to the dynamics of 
interacting elastic vortex lines in a disordered medium 
is still an open problem. The technical challenge con¬ 
sists in the accurate statistical summation of individual 
vortex-vortex interactions and pinning forces^which is 
akin to a many-body problem with long-range interac- 
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tions. For example, the critical currents may be es¬ 
timated based on energy and force balance considera¬ 
tions, see, e.g., reviews 8-10 only for the most simple 
cases, such as uniformly distributed weak pinning cen¬ 
ters. These estimates, however, do not describe the more 
practically important cases of multiple large-size defects 
occupying a considerable fraction of the material, as in 
most commercial high temperature superconductors. In 
such cases, the vortex pinning effects of different pinning 
sites are not simply additive, especially for high concen¬ 
tration of defects and in the presence of strong magnetic 
fields. In general, critical currents for such complicated 
defect structures can only be computed using large-scale 
numerical modeling. 

The rapid progress of high performance computing 
combined with sophisticated numerical methods and ma¬ 
terials synthesis provides a strong impetus to realize 
approaches with predictive power. A new area driven 
by this powerful combination is the materials genome 
initiative^that aims at predicting materials with de¬ 
sired properties: materials by design. An analogous 
eritieal-eurrent-hy-design concept can be envisioned in 
applied superconductivity, with the use of large-scale nu¬ 
merical simulations for the rational design of defect mi¬ 
crostructures in order to optimize the current-carrying 
capacity of superconducting materials for targeted appli¬ 
cations. The building blocks for this design are materials 
defects, which range from atomic to mesoscopic in size 
and their interaction with superconducting vortices. 

2. MODEL 

A well-established description of superconducting 
properties is based on the phenomenological time- 
dependent Ginzburg-Landau (TDGL) model.In this 
model, the superconducting state is represented by a 
complex-valued superconducting order parameter de¬ 
pending on time and coordinate, '0(r, t). Vortices enter 
as topological defects associated with a phase-winding 
of 27r in the order parameter field. This model auto¬ 
matically accounts for the interactions between vortices, 
between vortices and pin sites, for vortex-line elasticity, 
vortex cutting and re-connection, in other words, all the 
interactions that are relevant for describing pinning. The 
TDGL model can be reduced to a non-linear partial dif¬ 
ferential equation for the order parameter, given here in 
dimensionless form 

{dt + ija)'ip = s{r)'ip - + (V - iA)‘^^lj -h C(i*, t). (I) 

We use the temperature-dependent coherence length ^ = 
^(T) as unit of length, which allows to account for the 
temperature dependence of physical quantities in a com¬ 
putationally efficient way. The micro and nanostructure 
of the superconductor is modeled through the spatial de¬ 
pendence of the critical temperature Tc(r) contained in 
the prefactor of the linear term, ^(r). This prefactor is 
5(r) = I in the superconductor and is negative, ^(r) < 0, 


in the non-superconducting defect. A is the vector po¬ 
tential describing the magnetic field as B = curl A, fi is 
the electrostatic potential, and (^(r, t) is the temperature- 
dependent Langevin noise. We employed the so-called 
infinite-A approximation within which the vector poten¬ 
tial is fixed by the applied magnetic field. This approxi¬ 
mation is suitable for type-II superconductors with high 
Ginzburg-Landau parameter (A/^ ^ 100 for the experi¬ 
mental system) at high magnetic fields, when the inter¬ 
vortex distance is much smaller than the London pene¬ 
tration depth A and small spatial variations of the mag¬ 
netic field can be neglected. The details of the numerical 
simulations are presented in Appendix A. 

For realistic pinning microstructures. Equation (I) can 
be solved only numerically. This task proves to be com¬ 
putationally very demanding, and thus far, mostly meso¬ 
scopic or two-dimensional systems have been analyzed 
with this method, see e.g. Refs. 15-18. Recently, we have 
developed a stable implicit iterative solver for TDGL 
equation^^ based on the Jacobi method. This solver run¬ 
ning on high-performance general-purpose graphic pro¬ 
cessor units allows to simulate three-dimensional sam¬ 
ples large enough that finite size and surface effects be¬ 
come negligible.Predictions for macroscopic quan¬ 
tities such as the critical current are now possible, 
thereby enabling the concept of critical current by de¬ 
sign. We demonstrate this concept on technologically im¬ 
portant rare earth barium copper oxide (REBa 2 Cu 307_<5 
or REBGO) coated conductors and validate the results 
on samples with clearly identifiable pinning effects due 
to linear correlated defects in the form of self-assembled 
barium zirconate (BZO) nanorods and irradiation tracks 
introduced by heavy-ion irradiation. 


3. RESULTS AND DISCUSSION 

We present results on two sets of commercially avail¬ 
able REBGO coated conductors, one containing self- 
assembled BZO nanorods, the other without. The sam¬ 
ples were grown on Hastelloy substrates using metal 
organic chemical vapor deposition (MOGVD). The in¬ 
plane texture of a MgO buffer layer located between 
the superconductor and the substrate is achieved with 
an ion beam assisted deposition (IBAD) process.It 
has recently been shown that the addition of BZO can 
lead to the formation of self-assembled BZO nanorods 
inside the REBGO matrix during the synthesis of the 
superconductor.The nanorods grow largely along 
the c-axis of REBGO with typical diameters of ^ 5-10 nm 
and lengths of several hundreds of nanometers that is 
mostly limited by planar precipitates and stacking faults. 
Their addition results in substantial increases of the crit¬ 
ical current density with values reaching « 8MA/cm^ 
at 30 K and in a field of 9T applied parallel to c-axis. 
At high temperatures (T > 60K), the linear shape of the 
BZO nanorods is responsible for the pronounced peak in 
the field-angle dependence of Jc for applied fields aligned 
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Figure 1: a. Transmission electron microscopy (TEM) image of the unirradiated sample with self-assembled nanorods aligned 
along the c-axis. b. High-resolution TEM image of BZO nanorods. The nanorods are oriented along the c-axis, and the dark 
striations are Moire-patterns resulting from the imaging setup, c. Heavy-ion damage tracks revealed by Z-contrast scanning 
TEM in a cross-section view along the [100] direction with atomic resolution. The tracks are inclined at ~ 45° from the 
c-axis. In panels a-c the color scheme was chosen to enhance the contrast, d. Angular dependence of the critical current Jc{0) 
before and after irradiation in REBCO coated conductor samples at a temperature of 77 K and in a magnetic held of H = 1 T. 
The sample was irradiated with 1.4 GeV Pb^®^ ions at 45° from the c-axis. In the unirradiated sample, the presence of BZO 
nanorods induces a large peak when the held is applied parallel to the c-axis of the REBCO layer, and a small peak is also 
present for the a5-plane due to intrinsic pinning, RE 2 O 3 nanoplates, and stacking faults. In the irradiated sample, the peak has 
shifted toward the direction of irradiation, e. Representative conhguration of vortices and defects in the simulation. A snapshot 
of the isosurfaces of the order parameter is shown in red. The vortices follow the direction of the magnetic held applied along 
the c-axis. Columnar defects and nanorods (displayed shorter than actually used in the simulations for illustration purposes) 
are shown in grey. See the videos of vortex motion in the sample with short nanorods (magnetic hied B applied along the 
c-axis and at 45° to the c-axis) and in the sample with nanorods and tilted columnar defects (B along the c-axis and at 45° 
to the c-axis). f. Angular dependence of the critical current from simulations, normalized to the depairing current. Reference 
simulation with sample containing nanorods along the c-axis is shown in black. The addition of columnar defects at 45° from 
the c-axis shifts the peak as shown by the green dotted line. The addition of plate-like defects in the a5-plane modihes the 
angular dependence as shown in red. 


with the c-axis, whereas at lower temperatures, strain 
helds associated with the BZO nanorods and/or nanopar¬ 
ticles induce almost isotropic Jc-enhancements.^^’^^ 

An alternative approach to introducing linear defects is 
realized by irradiation with high-energy heavy ions.^^ For 
certain ranges of parameters — such as weight and energy 
of the ions, thermal and electric conductivity of the target 
material, and rate of energy transfer, continuous amor¬ 
phous tracks are formed in the target material. Typical 
track diameters fall in the range of 5-I0nm.^^“^^ These 
heavy-ion induced columnar defects are known to be the 


most effective vortex pinning centers for applied fields 
aligned with the tracks due to their ideal size for confining 
vortices, and their pinning effects have been extensively 
studied.In contrast to self-assembled BZO nanorods 
that are chemically inserted, irradiation induced colum¬ 
nar defects can be introduced at arbitrary angles and over 
a wide range of areal density,allowing for a conve¬ 
nient means of investigating the constructive or destruc¬ 
tive pinning action of several types of coexisting pinning 
sites. 

In this work, we present (i) a study of vortex pinning 
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in a REBCO coated conductor with self-assembled BZO 
nanorods before and after heavy-ion irradiation at 45° 
and (ii) a study of the cumulative effect of splayed heavy 
ion irradiation at angles of +45° and —45° to the c-axis in 
a REBCO coated conductor without nanorods. In both 
studies, we compare experimental and numerical results 
and demonstrate the validity of our numerical TDGL 
approach. Subsequently we explore a wider parameter 
range with numerical simulations to determine optimal 
pinning configurations. 


3.1. Sample with nanorods irradiated at an angle 

Eigure 1 summarizes our results on a REBCO 
coated conductor containing BZO nanorods into which 
columnar-shaped tracks have been introduced at an an¬ 
gle close to 45° from the c-axis using 1.4 GeV ^^^Pb ion 
irradiation. The unirradiated sample (Eigures la and 
lb) contains BZO nanorods with estimated density cor¬ 
responding to matching field ^ IT (meaning the 

effective defect density is equivalent to the vortex den¬ 
sity at a field of 1 T). After irradiation (Eigure Ic) sample 
contains 9.7 x 10^^ tracks/cm^ corresponding to matching 
field of B^ c ~ 2 T. (We distinguish the matching field for 
irradiated columnar defects, B^ from that of nanorods, 
N, and from for the total density of splayed columnar 
defects irradiated in two directions, These irradia¬ 

tion induced columnar defects have an amorphous core of 
^ 8 nm in diameter, similar to Pb-ion irradiation damage 
in YBGO single crystals and extend through the entire 
thickness of the film. Eigure Id shows a comparison of the 
angular dependence of Jc from transport measurements 
in unirradiated (black) and irradiated at 45° (red) sam¬ 
ples with self-assembled BZO nanorods. In these mea¬ 
surements, the current flows along the a+plane and is 
always perpendicular to the applied magnetic field. The 
field is rotated from the c-axis {0 = 0°) to the a+plane 
{0 = 90°). The irradiation-induced columnar tracks are 
oriented in such a way that they are perpendicular to the 
current flow. 

The angular dependence of the critical current Jc{0) 
conveniently characterizes the nature of the dominant 
correlated pinning centers. Pinning contributions from 
the inherent layered structure of the material, from in¬ 
plane nanoplates and from stacking faults create two 
strong peaks in for fields oriented along the a6-plane, 
at ±90°. As expected, the sample with BZO nanorods 
presents a pronounced peak in Jc near 
the expense of the in-plane critical current Jc{B\\ab).^^'^^ 
As a result, at 77 K and B = IT, the BZO nanorods actu¬ 
ally reverse the critical-current anisotropy from 2 to 0.5, 
defined as Jc{B\\ab)/Jc{B\\c) as can be seen by compar¬ 
ing the data labeled ‘pristine’ in Eigure 3d and the data 
labeled ‘nanorods’ in Eigure Id. 

Remarkably, in the sample containing both BZO 
nanorods and irradiation-induced columnar defects, we 
observe only a single peak in the critical current at 


0 ^ 45°, instead of the superposition of a Jc-peak due to 
the BZO nanorods at 6> 0° and a Jc-peak due to the 

columnar defects at ^ ~ 45°. Namely, the vortex pinning 
due to the BZO nanorods has been strongly suppressed 
by the introduction of the irradiation tracks. Such an ef¬ 
fect is a clear example of the non-additive nature of two 
types of pinning centers. 

When describing this situation theoretically, we mainly 
focus on the two types of pinning centers that domi¬ 
nate the angular dependence of the critical current and 
of which we have quantitative knowledge: nanorods 
and irradiation induced columnar defects. Eurthermore, 
we model intrinsic pinning by means of the Lawrence- 
Doniach model"^^ and a modified Laplacian describing 
anisotropy. In addition, we consider contributions from 
in-plane plate-like defects. The model of the tape nanos¬ 
tructure, that includes both the BZO nanorods aligned 
along the c-axis and the columnar defects inclined at 
0 = 45°, is presented in Eigure le. 

Our simulation allows for an unprecedented insight 
into the three-dimensional vortex dynamics as shown in 
Eigures le and 3e. We find that for large enough concen¬ 
trations of irradiation-induced columnar defects, vortices 
aligned with the BZO nanorods can easily jump from one 
nanorod to the next by sliding along the oblique colum¬ 
nar defects. This effect is directly responsible for the 
reduced critical current for magnetic field along the BZO 
nanorod direction. The resulting simulated angular de¬ 
pendences of Jc are shown in Eigure If and are in good 
qualitative agreement with the experiment. In particu¬ 
lar, they depict the non-additive effect of the two types of 
pinning centers and that the pinning action of the BZO 
nanorods is no longer visible. We therefore attribute the 
loss of pinning by the nanorods at 6> = 0° observed in the 
experiment to a ‘vortex sliding’ effect between nanorods 
induced by the oblique continuous columnar tracks, see 
insets in Eigure 2. In other words, the columnar defects 
tilted at 45° effectively reduce the nucleation energy of 
local vortex kink formation and thus, help vortices to 
slide from one nanorod to the next one. 

We typically observe that dominant extended defects 
define the critical current peak of Jc{0) at their tilt an¬ 
gle, while all other defects only slightly modify the shape 
of Jc{0), without adding additional peaks. Eor exam¬ 
ple, coated conductor samples contain a fair amount of 
stacking faults and a+plane precipitates that disrupt the 
columnar defects and the BZO nanorods. To elucidate 
their role, we simulate the system with and without in¬ 
plane, plate-like defects. The predominant effect of these 
planar defects is an enhancement of the critical current 
for in-plane magnetic fields, 0 = 90°. Comparing the 
dashed green and the red lines in Eigure If, one sees 
that the in-plane defects slightly reduce the peak in fields 
aligned with the columnar defects, 0 ~ 45°. 

Additionally, to obtain predictive results, we numer¬ 
ically analyze a wide range of concentrations of BZO 
nanorods and columnar defects. In Eigure 2 we show 
density plots of the simulated critical current as a func- 
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Figure 2 : a. Angular dependence of the critical current (normalized to the depairing current Jdp) for different concentrations 
of nanorods in a magnetic field of 0.02Hc2- One can see the prominent peak along the c-axis (0 = 0°) for a wide range of 
nanorod concentration (expressed in matching field ranging from ~ 0.025i7c2 to 0 . 1 i 7 c 2 . This corresponds to 

the volume fraction 5-20% occupied by nanorods inclusions. The peak critical current at ^ = 0° and ~ 0.07i/c2 is 

Jc = 0.14 Jdp. The white dotted curve shows the optimal concentration of nanorods for maximum critical current at a given 
angle of applied field. Inset: Sketch of the vortex pinned by a nanorod for B\\c. b. Angular dependence of the critical current 
for different concentration of columnar defects B^^c tilted at ^ = 45° and at a fixed concentration of nanorods (with matching 
held ~ 0.02iJc2 shown by dashed line in panel a). With increasing concentration of columnar defects, the peak shifts from 
^ = 0° to 45°. The near-optimal critical current corresponds to the angular range 0 = 30-40° and concentration of columnar 
defects B^^c ~ 0.02-0.04iJc2 (or volume fraction 6-11% occupied by them). Inset: Vortex sliding on tilted columnar defects. 


tion of angle 0 and concentration of the defects charac¬ 
terized by the matching field B^. Figure 2a shows that 
with increasing nanorod concentration (without colum¬ 
nar defects) a prominent maximum in Jc{0) emerges 
around 0 = 0°. For nanorod concentrations ranging from 
~ 0.025iJc2 to ^ 0.1iJc2 the height and angular 
width of this maximum are predicted to be almost inde¬ 
pendent of concentration. With further increase of the 
nanorod concentration, the Jc peak at J = 0° decreases. 

Figure 2a also shows that the two Jc peaks at J = ±90° 
attributed to intrinsic pinning from the material’s layered 
structure are suppressed when the nanorod defect den¬ 
sity exceeds B^ ^ « 0.04iJc2- The critical current for 
B\\c approaches its maximum at ~ 0.07iJc25 which 
corresponds to a volume fraction of 14% occupied by 
nanorods. We find that the mechanism of this suppres¬ 
sion is the same as described above for vortices jumping 
between nanorods by ‘sliding’ through the columnar de¬ 
fects, but this time it is the nanorods that create ‘bridges’ 
on which vortices can move from one intrinsic layer to 
the next. Experimentally, we note that the critical cur¬ 
rent for B\\ab due to intrinsic pinning above an over-all 
isotropic baseline level indeed decreases with increasing 
concentration of BZO nanorods, in agreement with our 
simulations. This trend holds even though the baseline 
Jc and additional pinning from planar precipitates — 
both not included in the current simulations — vary sig¬ 
nificantly between samples. 

Figure 2b shows the effect of increasing the concentra¬ 
tion of irradiation induced columnar defects for a fixed 
concentration of BZO nanorods of B^ ^ ^ 0.02iJc2- The 


intrinsic pinning peaks at ±90° vanish above a dose 
matching field of B^ c ~ 0.06iJc2- A similar reduction of 
intrinsic pinning was observed experimentally in YBCO 
films that contain inclined columnar defects.Further¬ 
more, the peak at ^ = 0° due to the BZO nanorods 
rapidly shifts to ~ 45° with increasing columnar defect 
concentration from B^ c ~ 0.02iJc2 to ~ 0.04iJc2- For 
even larger columnar defect concentrations, this latter 
peak becomes smaller as the effective barrier for vortex 
hopping between neighboring columnar defects decreases. 
This shift in the Jc peak when BZO nanorods compete 
with columnar defects oriented at a different angle clearly 
highlights the complex vortex dynamics in a mixed pin¬ 
ning environment. Although newly added defects can pin 
vortices, they may also facilitate vortex hopping or slid¬ 
ing between neighboring preexisting defects under spe¬ 
cific conditions. Hence the effect of these defects on vor¬ 
tices can be tailored by controlling the concentrations 
and mutual orientations of the defects. 


3.2. Pristine sample irradiated at two angles 

Non-additive effects can also occur between two iden¬ 
tical types of pinning centers. For instance, irradia¬ 
tion induced columnar defects oriented at ±45°, also 
known as splayed irradiation, may offer a route towards 
a marked reduction of the anisotropy of Jc at an over¬ 
all high value in single crystals and thin films.The 
REBCO sample without BZO nanorods was irradiated 
with 1.4 GeV ^^^Pb ions at both —45° and ±45° from 
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Figure 3: a. Z-contrast scanning TEM image along the [100] direction of the = 2T, ±45° irradiated REBCO tape. 

The damage tracks are at ±45° from the c-axis and extend through the entire thickness of the REBCO film. The pre-existing 
stacking faults and nanoplates are indicated by horizontal and vertical arrows respectively, b. TEM image along the ±45° 
irradiation direction showing a cross-section with the tracks appearing as black dots. c. High magnihcation image showing 
the cross-section of a damage track, 8 nm wide. The featureless dark region indicates an amorphous structure. The elliptic 
shape and the shadow are due to a small misalignment between the irradiation direction and the sample-cutting direction. In 
panels a-c the color scheme was chosen to enhance the contrast, d. Angular dependence of the critical current in pristine and 
irradiated (IT in each direction) REBCO tape samples at 77K. The critical current performance is significantly improved by 
the incorporation of splayed columnar defects, e. Isosurfaces of the order parameter in a model with splayed columnar defects 
and an external magnetic held applied along the c-axis, also see the videos for magnetic hied applied along the c-axis and at 45° 
to the c-axis. f. Numerically simulated angular dependence of the critical current for several densities R<j>,x of splayed columnar 
defects. 


the crystalline c-axis with a total dose of = 2T 

= IT in each direction). As expected, the Pb- 
ion irradiation generates perpendicularly crossed dam¬ 
age tracks at ±45° from the c-axis. The morphologies 
of the crossed columnar defects are presented in the Z- 
contrast scanning transmission electron microscopy im¬ 
ages of Figures 3a-3c. Similar to the REBCO hlms re¬ 
ported earlier,these hlms also contain a large number 
of correlated nanoplates and stacking faults parallel to 
the him plane. We apply the current along the a5-plane 
perpendicular to the splay and applied magnetic held. 
The changes induced by the crossed columnar defects on 
the angular dependence of Jc are shown in Figure 3d. We 
hnd that Jc is enhanced over a wide angular range. For 
the irradiated sample, the increase of Jc reaches 27% for 
B\\ab and 71% for B\\c. As a result, the critical-current 
anisotropy is markedly reduced down to 1.17 when de- 
hned as Jc{B\\ab)/Jc{B\\c), or 1.3 when more properly 


dehned as max^j Jc(6>)}/ min^j Jc(6>)}. 

We carried out TDGL simulation of the vortex dy¬ 
namics in REBCO containing splayed columnar defects. 
Figure 3e shows a representative snapshot of the arrange¬ 
ment of vortices and defects. The results of the sim¬ 
ulation for different concentrations of crossed columnar 
defects are shown in Figure 3f. The simulation quali¬ 
tatively reproduces the reduction of the anisotropy and 
the overall increase of the critical current observed exper¬ 
imentally in Figure 3d. At low concentration of splayed 
defects ~ 0.008F7c2)5 the simulation of the angular 

dependence of Jc reveals two maxima located at ±45° as 
might be expected for a simple addition of pinning effects. 
However, with increasing concentration the two maxima 
merge to form a single flat maximum around 0 = 0°, as 
in the experiment. Thus our simulations suggest that a 
fairly isotropic Jc can be achieved for relatively low de¬ 
fect concentration with a rather modest increase in Jc, 
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Figure 4: Angular dependence of the critical current for 
different concentrations of splayed columnar defects at ±45°. 
Similar to Figure 2a, a peak Jc = O.llJdp occurs at ^ = 0° 
for = 0.035i7c2, which corresponds to a volume fraction 
of 9.3% occupied by columnar inclusions. The white dotted 
curve corresponds to the optimal at given 6. 

whereas the maximum enhancement of Jc{B\\c) occurs 
at significantly higher defect concentration accompanied 
by a sizable reverse anisotropy, i.e. Jc{B\\c) > Jc{B\\ab). 

Figure 4 shows the simulated plot of Jc for different 
concentrations of splayed defects and angles of applied 
field. Again, as in Figure 2, one can observe two peaks 
at ^ = ±90° due to intrinsic pinning. The combination 
of columnar defects tilted at ±45° results in a broader 
peak near 0 = 0° compared to the peak in the angular 
dependence of Figure 2a for nanorods. The maximum of 
the critical current occurs at a matching field of B^ ^ = 
0.035iJc2, corresponding to a volume fraction of 9.3%. 

The optimal concentration of nanorods predicted for a 
REBCO conductor (Figure 2a) is approximately three to 
four times larger than the values observed in our experi¬ 
ments on currently available industrial samples contain¬ 
ing 7.5 mol-% Zr. Extensive recent work on research- 
scale samples has shown that the critical current den¬ 
sity increases strongly with Zr-content, that is, with the 
concentration of nanorods,and samples with 
Zr-doping up to 25 mol-% have been synthesized. Our 
simulations (Eigure 2a) are aligned with these studies 
and suggest that there is still room for further enhance¬ 
ment. In contrast, the predicted optimal concentra¬ 
tion of irradiation-induced splayed tracks in a pristine 
sample (Eigure 4) corresponds closely to the concentra¬ 
tion used in our experiments, and no further improve¬ 
ment is expected. However, we expect that the addi¬ 
tion of splayed columnar defects will decrease the critical- 
current anisotropy of the sample. 

Our large-scale TDGL simulations describing the joint 
action of different dominant pinning centers provide valu¬ 
able insights for the design of effective mixed pinning 
landscapes by elucidating the vortex dynamics responsi¬ 
ble for the enhancement of Jc in commercial high tem¬ 


perature superconductors. Visualizing the real dynamics 
is difficult, as few, if any, experimental techniques can re¬ 
solve both the microscopic magnetic structure and short 
timescale. It is however key to achieving higher and more 
isotropic critical currents in these high temperature su¬ 
perconductors. 

4. CONCLUSIONS 

In summary, we introduced the new critical-current- 
by-design paradigm. This paradigm aims at predicting 
the optimal pinning landscapes for maximum critical cur¬ 
rent in targeted high-temperature superconducting ap¬ 
plications. We illustrated this concept on technologically 
important rare earth barium copper oxide coated con¬ 
ductors and validated the results on samples with clearly 
identifiable pinning effects due to linear correlated defects 
in the form of self-assembled barium zirconate nanorods 
and irradiation tracks introduced by heavy-ion irradia¬ 
tion. The TDGL simulations elucidate the vortex dy¬ 
namics responsible for the non-additive behavior of the 
vortex-defect interaction in the presence of different types 
of correlated defects with various spatial orientations. In 
particular, we observed a vortex-sliding scenario, directly 
demonstrating the dynamics that can result from syner¬ 
gies between different types of pinning. Using this ap¬ 
proach, we also can predict the optimal concentrations of 
these defects for maximal critical current. As character¬ 
ization and models of realistic pinning centers improve, 
and as computational performance increases, the quan¬ 
titative prediction of a superconductor’s critical current 
from its microstructure now seems within reach. 
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Appendix A: Numerical simulations 
1 . Geometry 

We solve the time-dependent Ginzbnrg-Landan 
(TDGL) equation for a complex order parame¬ 
ter in a rectangular superconductor of finite size 
Lx:XLyXLc = 64^ X 256^ x 256^ = 218 x 870 x 870 nm^ 
with an underlying periodic numerical grid of size 
Nx X Ny X Nc = 128 X 512 x 768, where f is the 
superconducting coherence length. We used a mass 
anisotropy factor of 7 = 5 for YBCO. To simulate the 
experimental angular dependence, we apply a uniform 
external magnetic field B = 0 . 02 i 7 c 2 and rotate it in the 
^c-plane, see inset in Figure Id. 

In order to minimize surface effects and reduce the ge¬ 
ometrical factor in the angular dependence of the crit¬ 
ical current Jc{0) we compared two approaches. In 
the first approach we placed the whole superconduct¬ 
ing volume in a cylindrical domain surrounded by a non¬ 
superconducting metal and applied an external current 
along the cylinder, see Figure 5a. This approach shows 
a uniform angular dependence of Jc{0) for an isotropic 
system (i.e. with anisotropy 7 = 1) with randomly dis¬ 
tributed spherical inclusions. Conversely, a rectangular 
system with open (no-current) boundary conditions gives 
10-30% variations in Jc{0) under the same conditions. 
Using periodic boundary conditions in all directions (as 
shown in Figure 5b) also gives an almost uniform Jc{0) 
for the used benchmark system. In that case it is im¬ 
portant to use a sufficiently large system such that finite 
size effects get minimized, i.e., artificial self-pinning by 
vortices ‘wrapping’ around the system. We choose this 
second approach because it has 20 % more superconduct¬ 
ing volume and thus allows acquiring more statistics of 
the vortex dynamics. 


2 . Grid precision 

The TDGL approach can adequately describe the 
vortex-vortex and vortex-defect interactions on scales 
comparable or larger than the coherence length There¬ 
fore our approach does not require a grid resolution much 
finer than To determine the required numerical grid 
precision, we performed a series of simulations on a (32^)^ 
sample and varied grid precision h = LjN (distance be¬ 
tween grid points) from 0.1 to 1.0. We found that all grids 
with h ^ 0.6 give the same I-V curves and, correspond¬ 
ingly, the same critical current. We use = hy = O.bf 
and he = 0.33^ in the present article. 


10 



Figure 5: Reduction of the geometrical factors affecting angular dependence of the critical current, a. Superconducting volume 
in a cylindrical domain and surrounded by a non-superconducting matrix, b. (Quasi)periodic boundary conditions for the order 
parameter. 


3. Intrinsic pinning of YBCO 


The intrinsic layered structure of YCBO can be incor¬ 
porated in the time-dependent Ginzburg-Landau model 
in two ways: (i) by spatial Tc(r) modulation of layers 
parallel to the a 6 -plane and (ii) by controlling the effec¬ 
tive interlayer distance. The first method is characterized 
by three parameters: thickness of the superconducting 
layer, thickness of layer with suppressed superconductiv¬ 
ity, and the Tc in the latter. The second method is based 
on the Lawrence-Doniach model and uses the numerical 
discreteness of the mesh as natural pinning landscape. In 
this model the effective interlayer coupling is controlled 
by the parameter where 7 is the mass anisotropy 

factor and he = Lc/Nc. 

We found that these methods typically give similar re¬ 
sults that can be mapped to each other for a sample with 
periodic boundary conditions and without additional in¬ 
clusions. This means, we can obtain similar angular de¬ 
pendencies Jc (0) of the critical current with a fixed ratio 
Jc(90°)/Jc(0°) by tuning the parameters of the before 
mentioned methods. However, in reality the intrinsic lay¬ 
ered structure has a much smaller scale than the charac¬ 
teristic size of other inclusions. This makes the Lawrence- 
Doniach model more convenient to describe the intrinsic 
layered structure of YBCO for large-scale simulations as 
it, in a sense, integrates out the small scale intrinsic pin¬ 
ning behavior into the parameter 

Based on a series of simulations of the superconduct¬ 
ing sample with given concentration of nanorods and 
mass anisotropy factor 7 = 5 we found that a Lawrence- 
Doniach parameter ^/ 7 /ic = 0-6 gives the same ratio 
Jc(90°)/Jc(0°) of the angular dependence as in the ex¬ 
periment (compare Figure If versus the experimental an¬ 
gular dependence in Figure Id). We keep this empirical 
parameter in all simulations. 


4. Modeling of inclusions 

The inclusions are taken into account by the spa¬ 
tial variation of the critical temperature Tc(r), namely 
Tc(r) = Tc,s > T in the superconducting matrix and 
Tc(r) = Tc,i < T inside inclusions. Tc-modulation is 
naturally embedded into TDGL model by the variations 
of the coefficient ^(r) of the linear term in the dimen¬ 
sionless TDGL Equation (1). In this equation we use the 
coherence length ^ at a given temperature as unit length 
so that the prefactor ^(r) = 5s = 1 describes the super¬ 
conductor at this temperature. The non-superconducting 
inclusions with Tc,i < T can be described by a negative 
prefactor, s{r) = Cj = (Tc,i — T)/(Tc,s — T) < 0. Test 
simulations show that the value of the critical current 
does not depend on Si for any Sj ^ —1. Therefore, we set 
5 i = — 1 inside inclusions. 

Note, that this approach defines normal-metal inclu¬ 
sions subject to the proximity effect. Although real de¬ 
fects are believed to be insulating in their central region, 
this still gives a qualitatively correct result as the super¬ 
conducting material in the vicinity of defects gets dis¬ 
torted as well and therefore suppresses the order param¬ 
eter gradually. 

The prefactor 5 (r) and the Lawrence-Doniach parame¬ 
ter are the material parameters of our model. Therefore, 
the TDGL model can predict the optimal shape and con¬ 
centration of the inclusions, as long as 5 (r) [or Tc(r)] 
describe the inclusion materials realistically. 


5. Sizes of inclusions 

We determined the nanorods’ diameter from TEM im¬ 
ages of the first sample before irradiation to be about 4^ 
(with ^ ^ 3.4 nm) and used nanorod dimensions 4^ x 4^ x 
128^ with an average spacing of 16^, which corresponds 
to a matching field of ^ ^ 0.015iLc2 0.75) 

and volume fraction occupied by nanorods /n ~ 0.030. 
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Figure 6: Example of a simulated I-V curve shown in linear scale in panel a and in double-logarithmic scale in panel b. The 
critical current is determined by the intersection with a finite electric field level Ec- At currents lower than those shown in 
panel b, the vortex depinning events are so rare that the calculated electric field becomes very noisy. 


The angular dependence Jc{0)^ shown by the black line 
in Figure If, does not depend on the nanorods size along 
the c-axis provided they are longer than ^ 47 ^. Based on 
STEM images of the irradiated samples, the irradiation 
induced columnar defects have approximately the same 
diameter of 4^. The angular dependence Jc{0) shown 
by the green dashed line in Figure If was obtained for 
columnar defects inclined at 45° with a matching field of 
~ 0.032i7c2 ~ 1.6) and a volume fraction 

occupied by the columnar defects of fc ~ 0.083. In order 
to simulate the plate-like defects around each columnar 
defects as observed by TEM we added flattened ellipsoids 
of size 12^ X 12^ X ^ around the tilted columnar defects. 
The volume fraction /p ^ 0.017 of such defects reduces 
the peak near 45° and increases the anisotropy at ±90° 
as shown by the red line in Figure If. 

The angular dependencies shown in Figure 3f were 
obtained for matching fields from x ~ 0.008i7c2 

{B^^^/B « 0.4 and volume fraction /x ~ 0.021) to 
0.032i7c2 {B^,^/B ^ 1.6 and /x 0.084). 

Note, that the simulations do not account for the do¬ 
main structure of YBCO tape, such as twin boundaries. 


6. Finite voltage criterion 

The critical current is determined from the I-V curves 
by a constant electric field criterion with a threshold 
field of Ec = 5 x 10“^F^0: where the electrical field unit 
Eo = PabJo = (3v^/2)pa6Jdp ~ 2.QpabJdp naturally re- 
suits from the TDGL formalism, see details in Ref. 19. 
I-V curves are obtained by incrementing the external 
current from zero, while measuring the time average of 
the voltage across the system at each step, see Figure 6. 
Since vortex dynamics in the presence of pinning is hys- 
teretic, the I-V curve can depend on the external current 
ramping protocol and the averaging time. To ensure the 
validity of our procedure, i.e., calculation a unique I-V 


curve, we compared three different ways of obtaining it: 
ramping the current up, ramping the current down, and 
fixed current equilibration. All three approaches result 
in the same Jc- 

For YBCO the in-plane resistivity pah = 5 x 10“^ Et cm 
near the critical temperature and the depairing current 
^dp ~ 44MA/cm^ criterion can be estimated resulting in 
an electrical field threshold value of E^ ~ 0.3V/cm. To 
simulate dissipation levels as low as the industry standard 
of 1 /iV/cm, which corresponds to a very small hoping 
rate of vortices, one would need to simulate the system 
for a very long time, which is currently beyond reach 
in terms of computational cost. Typically, an experi¬ 
mental value for the voltage or resistivity is measured 
during 1 s. In the simulations we use 10^-10^ simula¬ 
tion time steps as averaging time, where each time step 
corresponds to 0.1 of the Ginzburg-Landau time and the 
Ginzburg-Landau time is of order 10“^^ s. This results in 
a reasonable numerical threshold value for which is 
still approximately five orders of magnitude higher than 
the industry standard, but since the I-V curve near the 
critical current has a sharp behavior (see Figure 6b), this 
does not translate into orders of magnitude differences in 
terms of Jq. In order to quantitatively characterize this 
threshold we roughly approximate simulated I-V curve 
near the critical current by a power law E ^ and esti¬ 
mate exponent n ^ 18, which means that the higher E^ 
criterion yields a critical current ^ 90% higher than the 
value one could expect by using the industry standard. 
Thus, our criterion may overestimate Jc by a factor of 
two. Such an agreement is remarkable in the field of vor¬ 
tex pinning. 

7. Langevin noise 

The simulation of realistic temperatures used in high- 
temperature superconductor experiments is a challenge. 
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The absolute temperature in the TDGL formalism is de¬ 
termined by the Langevin noise term in Equation (1). In 
Ref. 19 the concrete relation between temperature and 
noise characteristics is derived in detail. In this work 
we are limited to uncorrelated noise, which results in a 
technical upper limit in the temperature, which we can 
simulate- determined by the numerical stability of the 
integration scheme. To overcome this limitation, time 
and space correlations of the noise term are required. 
However, this is beyond the scope of this work and re¬ 
quires separate investigations. One can expect that a 
higher temperature leads to systematically lower critical 
currents, such that the qualitative behavior of our results 
is still valid (at experimental temperatures used here, the 
system is still not in the vortex liquid phase). 

Appendix B: Samples and measurement methods 

The two types of YBCO coated conductors used in 
this study (with and without BZO nanorods) were both 
grown using MOCVD^^ and have a critical temperature 
Tc ~ 92 K, coherence length ^(77 K) « 3.4 nm, London 
penetration depth A(77K) « 260 nm, and upper criti¬ 
cal field i7c2(77K) 28 T. The critical current densi¬ 

ties at 77 K and OT of the two types of YBCO films 
are ~ 410 A/cm-width, corresponding to Jc(77K) « 


4.1MA/cm^, while the depairing current density is es¬ 
timated to be Jdp(77K) 44MA/cm^. The critical 
current measurements were performed on four-contact 
bridges with a width of 85 /im, a thickness of 1.2 //m, and 
a length of 4 mm, patterned using optical lithography. 
The voltage-current characteristic was measured using 
a four-probe pulsed I-V method. The duration of the 
current pulse was 25 ms to minimize self-heating effects, 
and Jc was determined using the standard criterion of 
1/iV/cm. The characterization of the angular depen¬ 
dence of Jc was carried out at 77 K using a ^He cryostat, 
and the external magnetic field of 1 T was generated by a 
triple-axis magnet system which can rotate the direction 
of the field from —100° to +100° with respect to the film 
normal. 


Appendix C: Heavy-ion irradiation process 

YBCO bridges were irradiated at the ATLAS facility 
at Argonne National Laboratory with 1.4 GeV ^^^Pb ions 
to dose matching fields of 2T, either at 45° or at ±45° 
from the normal to the film. By using a low enough Pb 
ion beam current, heating effects are minimal and the 
transition temperatures of all the YBCO samples remain 
nearly the same after irradiation. 


